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ABSTRACT 

Non-monotonic features of rotation curves, and also the related gravitational effects 
typical of thin disks - like backward-reaction or amplification of rotation by negative 
surface density gradients - which are characteristic imprints of disk-like mass distri- 
butions, are discussed in the axisymmetric thin disk model. The influence of the data 
cutoff in rotational velocity measurements on the determination of the mass distribu- 
tion in flattened galaxies is studied. 

It has also been found that the baryonic matter distribution in the spiral galaxy 
NGC 5475, obtained in the axisymmetric thin disk approximation, accounts for the 
rotation curve of the galaxy. To obtain these results, the iteration method developed 
recently by the authors has been applied. 

Key words: gravitation: disk model, galaxies: kinematics and dynamics, individual: 
NGC 5475, mass function 



1 INTRODUCTION 

The analysis of rotation curves provides the most reliable 
means for ascertaining at least the gross distribution of grav- 
itati ng matter with in spiral galaxies, as Alar Toomre pointed 
out ljToomrelll963l l. In the same paper Toomre formulated 
a complete mathematical model of an axisymmetric and in- 
finitely thin disk rotating under its own gravity. The model 
offers a tool for determining the equilibrium mass distribu- 
tion directly from the rotation law of a highly flattened sys- 
tem, such as a spiral galaxy. It is assumed in the framework 
of the model, that orbits of stars, gas, etc., are circular, and 
that gravitational forces are balanced solely by centrifugal 
forces; the effect of pressure is thus ign ored. A few years ear - 
lier, following the work reported by tB urbidge at aLlflOSgf ). 
Brandt discussed less general situation of a very flattened 
system regarded as an assei nbly of oscu lating homoeoids 
compressed to a circular disk l|Brandtll960l l. All thes e papers 
follow ed many other pioneering ones referred to in ijBrandtl 
Il96(]h . 

The customary parametric few-component models re- 
late the mass distribution of baryonic matter mainly to lu- 
minosity measurements. The obtained amount of luminous 
mass is usually insufficient to account for the observed ro- 
tation of galaxies and a massive spherical dark halo is in- 
troduced as a remedy for the missing mass. The thin disk 
model with the surface mass density reconstructed mainly 
from rotation curves performs very well with strongly non- 
monotonic rotation curves, whereas the customary models 
have difficulties in explaining them. The thin disk model can 
easily account for high local gradients of rotational velocity, 
typically giving a lower total mass. In an extreme example. 



the rotational velocity of an outer galactic region, falling off 
faster than Keplerian, cannot be explained by the presence 
of a massive spherically symmetric halo. However, such a 
feature may be explained using a ffattened mass distribu- 
tion with a suitably changing mass proflle. 

The major difficulty in the practical use of the disk 
model lies with unambiguous mass density reconstruction. 
Unlike for spherical symmetry, the density depends on the 
assumed extrapolation of the rotation curve beyond the last 
measurement point (referred to in this paper as the 'cut- 
off radius'). This is a consequence of the nonlocal relation 
between the surface density and the rotational velocity of 
a disk-like system. For a reliable reconstruction of the disk 
surface density, the rotation curve has to be known globally 
(which, for observational reasons, is impossible) or at least 
out to its Keplerian falloffQ Therefore, the velocity mea- 
surements must be supplemented with independent data in 
order to constrain the sought mass distribution. This may 
be achieved, for example, by taking into account the amount 
of hydrogen measured in the out ermost galactic regions as 
was done in ijjalocha at abllioOSl l. 

Mass estimates of spiral galaxies are strongly model- 
dependent. Rotational velocities, when interpreted in a 
spherical dark halo model, may greatly overestimate galac- 
tic masses as compared to the disk model. It is thus im- 
portant to consider a phenomenologically acceptable model 
that gives a lower limit to the mass. The thin disk model 
may serve as such a reference model. 



1 ISofue fc Rubiiil |20o3) report on galaxies with Keplerian rota- 
tion curves at large radii 
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NGC 4736 is an example of a spiral galaxy with a ro- 
tation curve of which outer parts cannot be satisfactorily 
reconstructed when a spherical halo is assumed. By a sim- 
ple examination one finds that in the circular orbit approx- 
imation, the rotation curve of this galaxy cannot be created 
by a spherical matter distribution. By applying the thin disk 
model, with no constraining assumptions about the mass-to- 
light ratio profile, one can find a mass profile in this galaxy 
that perfectly conforms with its rotation curve, and agrees 
with the amount of hydrogen observed at large radii, beyond 
the cutofi^ radius ij Jalocha at al.|[^008l l. Only insignificant (if 
any) amount of dark matter is required, whereas the custom- 
ary models predict the galaxy to be dark matter dominated. 
In this paper we report also on another spiral galaxy, NGC 
5457, in which the amount of baryonic matter found in the 
disk model accounts for this galaxy rotation. 



1.1 The global thin disk model as a means of 
relating the rotation law to a mass 
distribution 

The main objection against the use of the axisymmetric thin 
disk as a model o f a spiral galaxy, is its instability with to o 
short time scale ijToomrd Il964l : l0striker~fc PeebledfToTi ). 
The simplest way to stabilize such a system is to add a suf- 
ficiently strong spherically symmetric potential. This is usu- 
ally done by introducing a massive spherical halo of dark 
matter, which is believed to surround spiral galaxies. How- 
ever, it should be noted that the galactic interior, irrespec- 
tively of its structure, produces at radii sufficiently large 
almost spherically symmetric potential that may stabilize 
the external circular orbits. 

Putting the stability issue aside, the disk model is still 
useful for approximate description of the gravitational field 
of a ffattened galaxy. Suppose that the rotation curve of such 
a galaxy represents the velocity of the streaming motion of 
matter rotating on roughly circular orbits in the galactic 
disk. The disk model associates with this rotation curve a 
formal surface mass density that may be considered as an ap- 
proximation of the column density of matter in this galaxy. 
To ensure its applicability, we use the model only for de- 
scribing galaxies with rotation curves breaking the condi- 
tion of spherical mass distribution at large distances from 
their centers. We assume that this non-sphericity can be at- 
tributed to flattening of these galaxies. In these regions we 
therefore expect the disk model to better reflect the oblate- 
ness of such galaxies than the models with massive spherical 
halo. What's more, the customary rotation curve modeling 
already assumes a resultant thin disk surface density as a su- 
perposition of the stellar disk component's surface density 
and of the column density of the spherical bulge compo- 
nent projected onto the thin disk's plane. Based on these 
premises the use of the single global thin disk for modeling 
the whole mass distribution in a flattened galaxy seems jus- 
tifled. The only new qualitative thing is that the effective 
disk becomes extended further out to the outer part in the 
case when the sphericity condition is broken there, that is, 
when the massive spherical halo cannot be introduced. 



2 THE SPHERICITY PROBLEM 

The velocity of test bodies moving on circular orbits in 
the equatorial plane of a spherically symmetric mass dis- 
tribution with a mass function M(r) is given by v{r) = 
•y/ GM{r)/r. By analogy with this relation we define the 
Keplerian mass function rv^{r)/G of a galaxy with a rota- 
tion curve v(r). For a galaxy with a massive spherical CDM 
halo, the Keplerian mass function, at least for large radii, 
should agree with the true mass function of the galaxy. 

Keplerian mass functions of some spiral galaxies decline 
or are non-monotonic in their outer parts, suggesting the 
presence of an extended, fiattened and massive subsystem 
in these galaxies rather than a spherical one. The spheri- 
cal dark halo may be excluded in favor of a large ffattened 
subsystem if the sphericity condition - 



for all r : {rv^irfj ^ 



(1) 



- is broken at larger radii (beyond the central bulge). In this 
case the global disk model seems more appropriate for re- 
construction of the mass distribution. Note, that the central 
spherical bulge may be represented by a column mass den- 
sity in the disk plane since the details of the internal mass 
distribution are not very important for the determination of 
the gravitational field in the distant regions, thus the use of 
the disk model in the central regions is also acceptable. 

Since there are many examples of disk-like mass configu- 
rations with their rotation laws satisfying the sphericity con- 
dition 111), violation of the condition at larger radii provides 
a strong argument for the presence of an extended fiattened 
component and against a significant spherical halo. For ex- 
ample, the Toomre thin disk of unit mass with the surface 
mass density <j{x) = (1 -I- a;^)"^''^/ (27r) (a; is a dimension- 
less radius in cylindrical coordinates), rotates differentially 
with velocity v{x) = x(l-\-x^^ ^^'^ ljToomrelll963l ). hence 

[xv'^{x)y — (l -I- s^) "'^^ > 0. As so, a spherically sym- 
metric mass distribution is possible with the same rotation 
curve. In the approximation of circular orbits, the spheric- 
ity condition is thus necessary, although not sufficient, for 
the presence of an extended and dominating spherical halo. 
Figure [U shows two rotation curves that break inequality JlJ 
at large radii. 



3 ROTATION IN THE GRAVITATIONAL 

FIELD OF AN AXISYMMETRIC FLAT DISK 

3.1 Basic equations 

Rotational velocity v{p) of circular motion of a test body 
around the symmetry axis and in the symmetry plane of the 
gravitational field of an axisymmetric thin disk is related to 
the disk's surface mass density cr(p) by 



v\pl 
4Gp 



V 



(2) 
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a) b) 

Figure 1. The upper plots show rotation curves of spiral galaxies a) NGC1808 and b) NGC4736. The lower plots show the corresponding 
Keplerian mass functions defined in the text. The unusual behavior of these functions (a mass function should be nondecreasing) suggests 
that matter distribution in these galaxies may be dominated by an extended and flattened galactic subsystem. 



see appendix lAI I The above integral exists in the 'princi- 
pal value' sense, hence the V symbol. The second part of 
the integral describes the backward interaction effect typi- 
cal for disks - the velocity of rotation on a circular orbit is 
influenced also by matter present outside the orbit. Gravi- 
tation of disks is therefore qualitatively different from that 
of spherically symmetric systems. Note, that for arbitrary 
cr > 0, integral ^ may attain negative values at some radii. 
For such radii circular orbits are impossible. 
The formal inverse of ^ reads 



1 



V 



P X 



X P - X 



■dX 



(3) 



This integral exists if lim pv{p) — 0, however in its deriva- 

tion more stringent assumptions were imposed on v{p) (see 
appendix [All. 



^ p is the radial variable in cylindrical coordinates and K and 
E are complete elliptic fu nctions of the first and second kind, 



llRvzhik fc GradsteinlfToSlh 



7r/2 



7r/2 



K(k) 



1 — siv? 4> 



E{k) 



In the approximation of cold, circular streaming mo- 
tion, relation ^ may be regarded as the definition of an 
effective surface mass density corresponding to the rotation 
curve of a fiattened disk-like galaxy. However, the answer 
to the question how accurately the resulting gravitational 
potential approximates the true one in this galaxy, would 
require separate studies, which are outside the scope of this 
paper. Nevertheless, it is still interesting to find out whether 
the effective surface density, that by construction accounts 
for the galaxy rotation, would also be consistent with other 
observational data available for the galaxy. Such a proce- 
dure led to coiisiste nt results for the galaxy NGC 4736 in 
ijjalocha at al.ll2008l l and for the galaxy NGC 5457 in this 
paper. 

Another relation between th e rotation velo city and the 
corresponding surface density is l|Toomrelll963 ) 



?(P) = ^ P 



dx P \p' 



+ /^.i^^Kldx 

dx X VX 



(4) 



where again V denotes a principal value integral. This for- 
mula became p opular mainly thanks to the handbook by 
iBinnev fc Trem aine (1987). Provided limp^o v'^ip) 0, for- 
mulas ^ and (|4]) are equivalent o-(p) = a{p) unless integra- 
tion gets cut off at a finite radius. However, the disadvantage 
in the practical use of the formula (jl} is that it contains the 
first derivative of v'^{p) which is subject to high observa- 
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tional errors that substantially enlarge the uncertainties in 
the determination of the surface density. 



3.1.1 An example 

To see how the basic equations of disk model work, we may 
apply them to the most familiar example provided by the 
Mestel disklfl Let a{p) = 2n'ap '^^^^ surface density. 

On substituting x ~ and x = ^~^Pi respectively, in the 
first and in the second integral in (l2|), we obtain 



M 



4G p 2-Kap 



{l+x)K(x) ~ E{x) _ ;^ 
x{l + x) ^ ^ 2 



that is, the rotational velocity is constant: v^{p) = 
Reversely, for a constant velocity Vo ~ , we obtain from 
((Sj, by similar substitutions as before. 



7r2Gp 



{l + x)K{x) - E{x) 
x{l + x) 



dx - 



M 



2tt Gp 2TTap 



The textbook equation (g} cannot be used for Mestel disk 
directly, since its rotation curve violates the requirement 
lim v{p) — (otherwise, one would conclude that a disk ro- 

fating with constant velocity has vanishing surface density). 
However, if we write v^{p) = v^Q {p — e) with small £ > 0, 
where 6 stands for the unit step function, then v'^'{p) = 



In 



VoS{p — e), and now we obtain from ^ (^{p) 

the limit e ^ we obtain the correct result a{p) = p 
for p > 0, since K(0) = ^. Again the formula ((Sj has proved 
more advantageous than the textbook formula (|4]|. 



3.2 Some consequences of the nonlocal 
correspondence between a and v 

Rotational velocities of disks may attain high values in the 
regions where the surface mass density decreases rapidly, 
and can be much larger than rotational velocities of spher- 
ically symmetric matter configurations of the same internal 
mass. To illustrate this, let us consider a sequence of unit 
mass disks with the surface mass densities 

a^{x)= ] + {2-{n + l)x-' 

+ (n- 6(1 -a:), n > 2. (5) 

O is the unit step function and x is a dimensionless radial 
variable. For n sufficiently large, cr„(a;)'s are almost constant 
and fall off rapidly to close to x = 1. Figure [2] illustrates 
the surface densities (t„(x) and the corresponding rotational 
velocities Vn{x). The functional sequence Vn{x) attains the 
theoretical limit 



y/^iKix)-Eix)), 0<X<1 
^■ix{K{x-^)-E{x-^)), x>l 



^ The example is simple but unfortunate since Mestel's disk is in- 
compatible with the mathematical consistency of the disk model, 
which was referred to in the appendix [aJ however, by doing cal- 
culations carefully, it can be still considered. 



which is unbounded at a; = 1. This signalizes that for 
large n and close to a; = 1, the rotational velocit y can be 
much larger than the Keplerian velocity vk {x) = M{x)/x, 
where M{x) is the mass function. This simple example il- 
lustrates that outside a highly oblate object, where densi- 
ties are already very small or negligible, the rotation curve 
may be still non-Keplerian even at distances from the object 
comparable with its size. 

Figure [3] shows another example of a differentially ro- 
tating disk. Its rotation curve resembles qualitatively rota- 
tion curves of some spiral galaxies. The first maximum in 
the rotation curve in figure Eljb) refiects the contribution 
from a central bulge, while the second maximum - from an 
exponential disk. The outermost data points on this curve 
would infiuence mainly the dark matter halo fit. If this model 
were to represent a real galaxy, the rotation curve could be 
measured only out to a; ^ 1, since beyond this radius only 
3 X 10~^ of the total mass would be present, thus amount 
hardly detectable. In effect, the total estimated mass pre- 
dicted by the customary models with dark halos, would be 
roughly equal to the Keplerian mass evaluated at the cutoff 
radius - the true mass in this example would be therefore 
overestimated by a factor of ~ 1.7. 

It is a well known fact, that external spherically sym- 
metric shells of matter do not infiuence the motions of stars 
on internal orbits, this is als o true when the shells of con- 
stant density are homoeoids dSinnev fc Tremainelll987l ). In 
general, however, the motions of stars are infiuenced by the 
presence of external masses. This backward reaction effect 
is particularly important for disk-like objects. To illustrate 
this, consider an axisymmetric system composed of a finite 
disk coplanar with an external massive ring. Figure [4] shows 
the rotational velocity of test bodies around the center of 
the system. In the presence of the external ring the veloc- 
ities are lower than those estimated based on the internal 
masses only. In the case when the external ring is sufficiently 
close or massive, such orbits may be impossible (the same 
is true if the internal mass is spherically symmetric). The 
presence of external fiattened shells of matter may help to 
explain the very declined parts of rotation curves of some 
galaxies. It is evident from figure |4] that the rotation of a 
very oblate system can change much with radius depending 
on the density gradients. In general, rotational velocities are 
amplified by negative and attenuated by positive gradients 
of the disk surface density. The gradients are not less impor- 
tant for the local value of the rotational velocity than the 
amount of matter. 

The elementary facts of disk gravitation discussed 
above, indicate that mass models of fiattened galaxies should 
be sufficiently general in order to take the disk specific effects 
into account. An important information about the mass dis- 
tribution encoded in the rotation curve may be simply over- 
looked by fitting a too much simplified parametric model, or 
by assuming the mass density in a luminous galactic subsys- 
tem to be proportional to its brightness profile. This in effect 
may lead to completely different qualitative and quantita- 
tive predictions for the same galaxy. A typical example of 
this situation is provided by the already mentioned galaxy 
NGC 4736. 
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Figure 2. The solid lines marked by a, b, and c, show respectively the elements n = 16, n = 50, and the limit n ^ oo of the sequence 
of unit mass disks with the surface mass densities cr„ defined in equation Jsj, and the corresponding rotational velocities. The rotational 
velocity in the equatorial plane of a ball with the same mass function as that of the limiting disk is marked by k. 
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a) b) c) 

Figure 3. a) The surface mass density of a unit mass disk: (j{x) = (48e-6400x2 _,_ 4g-400:r2/3 ^ g-iooa^^/s) + 

2^ {2 — YJx^^ + 15a;^^) 0(1 — x). b) The corresponding rotation curve (thick line) and its Keplerian asymptote (thin line), c) The 
true mass function (solid line) and the Keplerian mass function calculated from the rotation curve (dotted line). Although there is only 
3 X 10~^ of the total mass outside the radius x = 1, the rotation curve is still non-Keplerian and the Keplerian mass is almost twice as 
big as the true mass. To explain this rotation curve at large radii, one would normally add a massive spherical halo by fitting it into the 
rising part with the help of the least square method. In this way additional 'missing' mass would be introduced to the system to account 
for its rotation, while in fact this mass is not present at all. 




02468 10 02468 10 02468 10 



a) b) c) 

Figure 4. The influence of a massive external ring on the value of rotational velocities of test bodies moving on the internal orbits, a) 
The surface mass density of an axisymmetric system composed of an internal disk coplanar with an external ring; b) the corresponding 
rotation curves (and their Keplerian asymptotes) in the situation with (2) and without (1) the external ring; c) the mass function. Note, 
that in the presence of the external ring the rotational velocities in the vicinity of the internal disk are lower than those in the situation 
without the ring. 
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4 THE CUTOFF ERROR IN THE DISK 
MODEL 

As follows from equation the mass function of a disk- like 
system is 



o-{p) = o-{p), the corresponding cut off integrals are differ- 
ent: (t_r(/9) 7^ aFt{p), thus also different are the cutoff errors, 
Aanip) / AanXp). 

A cutoff error of mass determination of the internal disk 
of radius R, and corresponding to Acr_R, reads 



M(p) = — I dp\r I v'{px)n(x)Ax 



AMr 



Rvl 2 



G 



— I xAx 



dg Au^{^)^Jx 



where 



x^ - 1 



The total mass contained within radius p can be therefore 
determined only when a rotation curve is known globally, 
whereas the rotation curve of a real galaxy is known only out 
to a finite distance from its center (we call the distance the 
cutoff radius). This is the mathematical reason why the sur- 
face mass density and the mass function of a flattened galaxy 
cannot be determined uniquely from the measurements of 
its rotation curve only. This situation is drastically different 
from that for a spherically symmetric system in which case 
the mass function at a given radius is related directly to the 
local value of rotational velocity: M{r) = rv^{r)/G. 

Let R denote the cutoff radius. The integration in ([3]) 
gets naturally cut off at p = R, and for p < R we may 
approximate ^ by 



+ / v'ix) 



P X 

E(^) 

X - P 



p<R. (6) 



Since o-{p) ^ a_R(p), expression ^ gives the surface den- 
sity underestimated with respect to JSj by a cutoff error 
A(jr{p) = (j{p)~ (tb{p): 



where x = p/R, u{x) — v{xR)/vr and vr = v{R). 

Let us consider a theoretical situation of a disk-like 
galaxy with a rotation curve known to be almost Keple- 
rian outside a cutoff radius R. One then expects all of the 
galaxy mass M to be contained in the disk of radius R. Then 
also M ^ G~^Rvr and u{x) « for a; > 1. By using 

crR{p) to estimate the mass distribution in this particular 
situation, we obtain for < x < 1 



47r 



^■\AaRip)l^m^il^Pl£M^l^, (10) 



thus the resulting total mass is greater by |AMfl| = (4/7r — 
1)M ~ 0.27M than the actual mass M. The cutoff error 
|Acrji(p)| in lfTO|l grows with the distance from the galactic 
center. 

This example leads to the following criterion: provided 
that most of the disk mas is contained inside a radius i?, 
the cutoff error is insignificant in the region where cfr^p) 
satisfies the inequality 



(^Rip) > 



(11) 



7r2 RG' 

Otherwise, the approximation ^ is not reliable, and this 
occurs at larger radii, close to R. In the same theoretical 
situation, the cutoff error corresponding to the integral ([GJ, 
is given for < j: < 1 by 



RG . , , 1 
-2- • ^^r{P) = 



dc ^(f) 

5 C^-x'^ 



(12) 



AaR{p) 



Gn'^ 



ix) 



2 2dx>0, p<R. (7) 
X^ - P^ 



This error is completely unknown. With the Toomre integral 
([U cut off at p = J? the situation is even worse. In this case 
we have for p < i? 



dx P \P' 



p<R, 



(8) 



thus even the sign of the resulting cutoff error AaR^p) = 
(j(p) — o-r{p) is unknown, where 

oo 



since it depends on a weighted slope of the unknown part 
of the rotation curve. It should be stressed that although 



hence 
1 

4-K 



E{x)ln{l- 



27r2a 



x^) RG 
— < 



AcrH(p) < - 



ln(l-a:^) 
47r x^ 



where the inequalities 1 = E{1) < E(x) < E{x/£) < E{0) = 
7r/2 for X € (0, 1) and ^ > 1 have been used. The rightmost 
expression is divergent as x ^ 1, but for x < Xc ~ 0.63 it 
is still less than 7r~^, thus A(7_r(p) and |Acr_R(p)| are compa- 
rable for p < 0.63-R. In general, however, nothing is known 
about the mass distribution beyond R, therefore the cutoff 
error can not be even roughly estimated - for if, e.g., outside 
a cutoff radius R there existed a ring of dark matter con- 
centric with the disk, like in figure [H then the error could 
be arbitrarily large depending on the undetectable annulus' 
mass. 

In figure [5] compared are different approximations of 
the surface mass density obtained with the help of various 
methods, for the rotation curve of a model disk with a known 
mass distribution. To illustrate the infiuence of the spatial 
extension of the rotation curve measurements on the accu- 
racy of the reconstruction of the mass distribution, in figure 
E] presented is a similar analysis for various samples of the 
same model rotation curve that was cut off at different radii. 
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a) b) c) 

Figure 5. a) The rotation curve (and its Keplerian asymptote) of a disk with the surface mass density a{p) = 
^ {^e~'^f + j2e~f + Z^'*^). b) The surface mass density reconstructed with the help of different methods from the disk's 

rotation curve cut off at i? = 30: o - the exact profile, sf - the one calculated from equation Jsj, r - from equation JgJ, x - the sf 
surface density corrected for the contribution from the Keplerian tail according to equation ||10[I . sm - the surface density obtained from 
a spectral representation of the rotation curve (the internal mass was underestimated by 12%) - the surface mass density was calculated 
in this case from equation ||B4|| . The horizontal line represents the critical value of the surface mass density defined in equation llllll . 
the abscissa of the vertical line is ~ 0.63-R, which is the radius, to which sf and r may be considered comparable, c) The corresponding 
mass functions. 



As follows from the resulting mass functions, the cutoff error 
can be relatively large even in a distant empty region, where 
the rotation curve is not yet Keplerian, that is, when the cut- 
off error can not be estimated yet and taken into account. 
For example, in the surface density diagram in figure [5] the 
line representing aR{p) corrected for the cutoff error lflO)l . 
does not overlap with the exact value of the surface density. 
Even though 99, 9% of the total mass is concentrated in the 
disk of radius R = 30, the rotation curve is still not Keple- 
rian outside this radius. At this radius equation ^ overes- 
timates the true mass by 33%, whereas another estimation 
based on the spectral approximation given in appendix [B] 
underestimates the true mass by 12% only. 



5 MASS-TO-LIGHT RATIO AND THE 
RESULTS FOR NGC 5457 

The mass distribution of the luminous matter in differ- 
ent galaxy subsystems is customarily assumed to be pro- 
portional to the corresponding local brightness via con- 
stant mass-to-light ratios (frequently the blue band is cho- 
sen). This arbitrary assumption ignores the dynamical in- 
formation about the mass distribution encoded in the non- 
monotonic features of rotation curves. This in turn, may 
result in wrong distribution of matter among various galac- 
tic components. Especially important are the distant regions 
where a massive spherical hallo is always assumed, whereas 
in some cases a disk component of baryonic matter may 
be still more pronounced. Th e best illustration o f this is 
again the galaxy NGC 4736 ijjalocha at al.ll2008l ). There- 
fore, it seems better not to correlate the M/L ratio with the 
mass distribution and derive it as a point-dependent func- 
tion from the surface density determined in another way. 
Then one would be able to trace the presence of different 
star populations as well as the dark matter distribution. 



NGC 5457 (Hubble type Sc I) is one of the largest 
nearb y spiral galaxies, thus also well studied |Bosma at al.l 
Il98lh . The behavior of the local M/L ratio for this galaxy 
inferred from the Ks band differs from that in the B band, 
c.f. figure[7llQEI The reconstruction of the surface mass den- 
sity strongly depends on the chosen luminosity band - in the 
most part of NGC 5457, the local M/L ratio calculated from 
the B band is significantly greater than that from the Ks 
band. What's more, low-mass, lower temperature stars make 
up the most of galactic disk. Therefore, the global M/L ratio 
in the infrared band, like I, K, i, etc, is the best indicator of 
the presence of dark matter in spiral galaxies, with M/L < 2 
being typical for normal stars. 

Although the galaxy NGC 5457 is very large, its rota- 
tion curve has been determined only out to 14 kpc (NGC 
5457 has an asymmetric Doppler image). The rotation curve 
breaks the sphericity condition HI), c.f. figure[7l thus, for the 
reasons discussed in the previous sections, we use the global 
disk model. To obtain the global mass distri bution in this 
galax y we applied the method presented in ijjalocha at al.l 
I2OO8I ) that combines through iterations the measured ro- 
tation curve together with gas measurements far from the 
galactic center. The method gives a self-consistent global 
mass distribution in the disk. Self-consistency means that a) 
the rotation curve calculated from the global surface mass 
density with the help of integral ^ overlaps with the mea- 
sured rotation curve, b) the calculated global surface mass 
density outside the cutoff radius overlaps with the amount of 
hydrogen and helium visible in the outer parts of this galaxy 

* These results were reported in part at the XLVII Craco w School 
of Theoretical Physics in Zakopane jjalocha at al.|[2"o07^ . 
^ The measurement data for NGC 5457 were t aken from: th e HI 
surface density fc the L b luminosity profile llBraunlll99'it. th e 
galaxy distance lljarretl), the H2 surface densitv I^Young^ll991^ . 
and the rotation curve I Sofue. www pag^ V 
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Figure 6. The analysis of disk mass distributions which are superpositions of: three gauss profiles discussed in figure [5] (upper plots), and 
of two exponential profiles (lower plots). The analysis illustrates the influence of the cutoff radius on the accuracy of reconstruction of the 
original mass distribution. The surface densities were calculated with the help of integral JSJ (the lines marked by primed letters) and with 
the help of the spectral representation ||B4|| of the rotation curves (the lines without primes), a) Various pieces of the corresponding exact 
rotation curves obtained by cutting off external parts at different radii (the curves are shown in normalized variables (x, u) explained in 
the text), b) The reconstructed surface mass densities and the original density (denoted by 'o'). c) The corresponding mass functions, 
and the exact asymptotically unit mass function. 



and c) the calculated global rotation curve and global sur- 
face density are mutual transforms of each other according 
to integrals (HJ and ©. 

The global characteristics of the analysis for NGC 5457 
are given in table [T] In combination with the observed lo- 
cal brightness A, the determined surface density a yields the 
corresponding local mass-to-light ratio profile cr/A. In the 
Ks band, cr/A increases near the disk edge, c.f. figure [T] 
A simila r increase was observed b y Sofue for several other 
galaxies l|Takamiva fc Sofudl200(]| ). This rise would seem at 
first as the signature of the presence of dark matter in the 
outer part of the galaxy. Indeed, if stars were to comprise 
the whole baryonic mass of the galaxy, the a/ A ratio would 
be always that of stars, even in the limit A ^ 0. This is not 
the case for NGC 5457, however. The increase of the local 
mass-to-light ratio in this galaxy is caused by the presence of 
large amounts of hydrogen beyond the stellar disk. If there 
is some gas, eg. hydrogen and helium, then the local mass- 
to-light ratio profile (altars -f crH,Hc)/A diverges as A ^ 0. 
But the gas contribution should be subtracted as it is nonlu- 



minous in the visible and infrared bands. As is seen in figure 
[2l after subtraction of the hydrogen and primordial helium 
contribution, the local mass-to-light ratio profile stops to 
increase at the disk edge. Also in the B band the ratio has 
decreased. This fact, together with the low mass-to-light ra- 
tio in the band, shows that the CDM halo is not needed 
to account for this galaxy's rotation. Furthermore, the fact 
of breaking the sphericity condition, practically excludes a 
massive and spherical halo of CDM in NGC 5457. 



6 CONCLUDING REMARKS 

For the accuracy of determination of the mass distribution 
and its features in the disk component, local non-monotonic 
properties of rotation curves are important. They are not 
less important than the value of the rotational velocity. In 
particular, a local increase in the velocity field may be caused 
simply by a local decrease of matter density in the disk com- 
ponent and not by the increase of density in the spherical 
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c) f) 

Fi gure 7. The results for spiral galaxy NGC 5457 (MlOl) obtained in the thin disk model with the use of the iterative method developed 
in jjalocha at aljEoOSi '). a) the rotation curve calculated from equation |[2j for a surface mass density marked by u in figure b); b) the 
global surface mass density found by applying the iteration method, the hydrogen surface densities for HI , H2, and the surface luminosity 
profiles (A) (K and B bands); c) the Keplerian mass function Hll|l for the observed part of rotation curve of NGC 5457; d) the mass-to- 
luminosity ratio profile (o"/A) in Kg filter for total surface density and with hydrogen HI, H2 and helium subtracted (since they are not 
seen in the Ks band); e) comparison of local mass-to-luminosity profiles (cr/A) in B and Ks bands (hydrogen and helium subtracted); 
f) the mass-to-luminosity ratio profile in B filter for the total surface density and the other one with hydrogen HI , H2 and helium 
subtracted 
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Mtot 


1.13 X lO"Af0 


Mhi 


1.5 X W'^°Mq 


Mh^ 


3.4 X IO^Mq 


Lks 


7.45 X IO^OLq 


Lb 


2.49 X 10^0 Lq 


Mtot/Lks 


1.52 


Mtot / Lb 


4.54 


Mstar/Lks 


1.27 


Mstar/Lb 


3.8 



Table 1. The other results for NGC 5457. Mj^^q = 3.45, Mstar = Mtot 



I [Mhi - 



' A/h, ), D = 7.2Mpc, inclination ISdeg 



component of a galaxy. The important information about 
mass distribution carried by rotation curves may be simply 
overlooked if one assumes that the local mass distribution is 
proportional to the local brightness. In some cases the dark 
matter halo is introduced simply because the amount of lu- 
minous matter with the constant mass-to-light ratio found 
by the least square fitting method, cannot account for the 
rotation at large radii. The example of th e spiral galaxy 
NGC 4736 discussed in | Jalocha at al.lbOOSi ') very well illus- 
trates this observation. This shows also that determination 
of mass distribution in spiral galaxies is model dependent. 
For more reliable predictions it is thus important to have a 
generic model, more fiexible and general than the paramet- 
ric models with smooth profiles and constant mass-to-light 
ratios. It seems that the global thin disk model provides such 
a generic model, at least for fiattened galaxies with rotation 
curves breaking at large radii the sphericity condition. 

Due to observational cutoff in rotation data, the sur- 
face mass density reconstruction in the galactic disk (and 
consequently in the whole galaxy), is subject to high un- 
certainties, which increase with the distance from the cen- 
ter. Only in particular situations of rotation curves with al- 
most Keplerian tails, the cutoff errors can be eliminated and 
taken into account. The observationally determined part of 
a given rotation curve may be explained by various inter- 
nal disk mass distributions depending on how the rotation 
law has been extrapolated beyond the cutoff radius. This 
uncertainty can be removed if additional constraints on the 
mass distribution in the outer parts, not covered with rota- 
tion measurements, are taken into account. In this respect 
the hydrogen distribution in the outer regions can be used. 
We have established for galaxies NGC4736 and NGC5457 
that the observed amount of baryonic matter accounts for 
rotation of these galaxies in the approximation of the global 
disk model. 
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APPENDIX A: BASIC EQUATIONS OF THE 
AXISYMMETRIC THIN DISK MODEL 

We assume an axisymmetric mass distribution over an in- 
finitely thin disk. The cylindrical coordinates (p, (j), z) are 
chosen such that the disk overlaps with the z — G plane. 
Gravitational potential $ is z-symmetric, hence it suffices 
to consider the half space z > 0. In this region the complete 
function space of axisymmetric, everywhere bounded and 
vanishing at infinity functions is spanned by the base solu- 
tions Jo{ijjp/L)e~'^'^^^ of the Laplace equation, with u > 
and L being a length scale. Linearity implies that the most 
general solution in this function space can be represented as 
a superposition 



-2^1 



exp 



-a;— I Alu, 



(Al) 
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VL is a velocity scale, and the spectral amplitude a is dimen- 
sionless. A dimensionless distance x and a velocity function 
u{x) are defined by 



u(x) — 



v{Lx) 

VL 



The discontinuity of at 2 = is interpreted as an in- 
finitely thin layer of mass spread over the disk's plane with 
the surface mass density 



hence 



2ttG 
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oo 



(A2) 



Note, that is the spectral amplitude oi(j{p) (the inverse 
Hankel transform) 



GL 



(t{ijj) — — — I xa{Lx)Jo{uJx)d2 

^ T 



(A3) 



The condition for coplanar, circular and concentric orbits in 
the z = plane reads v'^(p) = p9p$(p, 0) for all p. Hence, 
having formally differentiated ijAip . one gets for z = Q 



u {x) = 2-KX J uja{uj)Ji{ojx)doj, 



and the resulting inverse relation 



a-{Lu) — — / u {x)Ji{ijjx)dx. 



(A4) 



(A5) 



On substituting IIASP to l|A2p . integrating by parts, provided 
that v{p) is differentiable, that u^(p) ^ as p — > 0, that 
w^(p)/V£ ^ as p — > oo|j and using the upper integral in 
footnotqj, we arrive at the Toomre integral 



^.iA'f^^ldx 
dx P \p' 



dx X VX 



(A6) 



The integral should be understood in the principal value 
sense (the symbol V indicates this kind of integration), thus 
the singularity of the elliptic function K{y) at y = 1 does 
not contribute. The equivalent form of the Toomre integral 
(with the same assumptions about u(p) as above) reads 



E(^) 



P X X P - X 



(A7) 



^ These assumptions assure at the same time the existence of 

integral |)A5|| and the suitable boundary conditions required for 

the equivalence of integrals 1IA6I1 and llATIl (flat rotation curves 

are incompatible with these conditions). To assure at the same 

time mutual invertibility of integrals llA4ll and llASIl one needs a 

stronger condition that lim ^v^{p) — > 0. 

p — '00 



J dujJo{u!x)Jo{uj^) ■■ 



00 

J dujJl{LtJ^)Jl{LtJx) ■■ 



cf. l|Rvzhik fc GradsteinlfTisih . 



IVX \ X / 



<x 
<x 



2 / ^ 

The integrand has an ordinary pole 2p(x-p) X. = P, which 
is easily integrable numerically and dominates the singu- 
larity of the elliptic function K. The above equation can be 
derived analogously like the Toomre integral, however it can 
be also proved by integration by parts the integral l|A6P (in 
order to use the relevant theorems of integral calculus, the 
integral l|A6P should be first rewritten as + J^^, then 
integrated by parts, and finally the limit £ ^ should be 
taken). 

The inverse of integral l|A7P reads 



AGp 



V 



-(x)^r^dx 



(AS) 



Indeed, by substituting l|A3ll in IIA4p one gets 



u^{x) 



2n 



Ji{ujx)duj I 



<^{Ly)dy (yJi(ujy))dy 



where the identity ui'^dy {yJiiujy)) = yJo{u)y) has been 
used. Next, on integrating by parts, with the assumption 
that p^a{p) — > as p — > and ^/p(y{p) ^ as p ^ 00, we 
obtain 



271- / ^d^d^^am) / Ji{Lj^)Ji{LJx)dLj 



u^{x) _ GL 



The rightmost integral is given in footnoti^. On integrating 
by parts again, with the same assumptions about the limit- 
ing behavior of a{p) as before, one arrives at the final result 
(IA8II . When doing these calculations, the same integration 
splitting and taking the limit, as in the derivation of ljA7p . 
applies. 

A simple analytic example using the above formulas and 
illustrating the equivalence of l|A6P and l|A7P was given in 
section 13.1.11 



APPENDIX B: A DISCRETE SPECTRAL 
REPRESENTATION 

For the observational reasons, a rotation curve «(p) of a 
galaxy is known only for radii p G (0, -R), R is the data 
cutoff radius. Let's define x = p/R, u{x) — v{Rx)/vr with 
VR = v{R). The function u'^{x)/x is thus defined on the unit 
interval x £ (0, 1) and m(1) — 1. Under some assumptions 
the function u'^{x)/x can be represented in this interval as 
a series of a complete set of orthogonal functions on this 
interval. Although we may choose any complete set of func- 
tions satisfying the appropriate boundary conditions on this 
interval, best suited for the disk symmetry are cylindrical 
functions. liu'^{x)/^ is integrable in the interval x G (0, 1) 
then i^(x)lx can be re presented in th is interval as a se- 
ries of Bessel functions IILebiedievI Il953l l. To conform with 
the more general integral representation l|A4p we choose the 
Bessel function of the first order, then 
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k 



GkJi{uJkx), 0<x<l, Jo(tJfc)=0. (Bl) 



The summation is taken over all positive zeros uj^ of the 
Bessel function of the zeroth order Jo (a::). On multiplying 
both sides by x Ji{uim x) and integrating over x G (0, 1) one 
gets 

1 

Ji(u;k) iOk J Ji(c^fc) 



The countable set of (t^'s suffices to describe «(p) completely 
only for p < R. However, this information is insufficient to 
obtain the continuous spectrum a{u}) of the true global o-(p), 
thus even for p < R one is unable to determine cr(p) from 
the observed part of rotation curve. This can be done only 
approximately. For example, assuming that v{p) = for 
p > R, we obtain from l|A5P 



1 _ LjJo{Lo)Ji{u;k) 

> /Ti. 



^2 -j^ ^ 



(B4) 



where a^'s have been calculated with the help of ljB2p . Let 
pi > be the smallest radius at which E(p) becomes nega- 
tive. By definition one sets in this approximation of the true 
(t(p) that cr(p) (Ja(p) where o-a(p) = for p > min(pi, R), 
and oa{p) = S(p) for < p < min(pi,ii). Then close to the 
center, the integral l|A8P with cr(p) — <ja{p), quite well ap- 
proximates the original rotation curve that was cut off at R 
(c.f. section [4] for more details). 



Since Joiu)) ~ Ji(ijJk)(i^k — <^) when u « oj^, the spectrum 
is continuous. Another extension of v(p) beyond R arises 
by assuming that expansion IjBip is formally valid also for 
p > R, then we obtain from l|A5P a discrete spectrum 



1 _ 5{lo — ujk) 



2-K ^ 

k 



(B3) 



It should be clear, that in a similar way, one can construct 
from the d-fc's set, infinitely many other spectra They 
will correspond to different global surface densities giving 
the same rotation law for p < R and different for p > 7?. In 
particular, this is the case for the global surface mass den- 
sities obtained by taking the inverse transform of the above 
two spectralfl However, the closer to the galactic center the 
less the two mass distributions differ from each other. In- 
deed, the error of determining o", which is calculated from 
equation ([9|) in the limit of small radii, reads 



limAa«(p)=^^ 



1 

2^ 



\0 



u{l) = 1, 



thus should be much smaller than the characteristic density 
scale v\/{RG) (as was shown in section |4] it is about 8% of 
the scale for Keplerian falloff), while the central density is 
usually much greater than this scale. 

With this o(ijo) we obtain from ljA2p the following series 



The fact that a function defined on the whole real axis, when 
restricted to a finite interval can be represented by many discrete 
and continuum spectra, should not astonish. To give a simple 
example, consider a function such that f(x) = sin x for x S (0, tt) 
and f{x) = elsewhere. Then 

27r(l — uj-') I 

are distinct spectral representations of this function on the inter- 
val (0, tt), that is, the inverse transform j f(u])e^'^^Auj con- 
verges to sin(x) in this interval for both the spectra. Note, how- 
ever, that outside the interval the integral converges to different 
functions. 



